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ABSTRACT 


The purpose of this effort has been to provide shortwave spectral radiance 
models for use in the spectral correction algorithms for the ERBE Scanner Instru- 
ment. The required data base was delivered to the ERBE Data Reduction Group in 
October, 1984. It consisted of two sets of data files: (1) the spectral bidirec- 
tional angular models and (2) the spectral flux models. The bidirectional models 
embody the angular characteristics of reflection by the Earth-atmosphere system 
and were derived from detailed radiance calculations using a finite difference 
model of the radiative transfer process. The spectral flux models were created 
through the use of a delta-Eddington model to economically simulate the effects 
of atmospheric variability. By combining these data sets, a wide range of radi- 
ances may be approximated for a number of scene types. 
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I. INTRODUCTION 


The algorithm being used to compute the spectral correction coefficients 
for the Earth Radiation Budget Experiment (ERBE) Scanner instruments requires 
a priori models of the spectral radiance at the top of the atmosphere (TOA) 
being sampled by the instruments. The purpose of this work therefore has been 
to provide those shortwave models (0.25 - 5.0 jam) through the use of well- 
established radiative transfer techniques combined with the most realistic 
input data available. The approach decided upon has been to use a detailed 
finite difference code to compute the TOA radiances to determine spectral 
angular bidirectional models for each of the surface and sun angle combinations 
with "mean" atmospheric conditions. At this point an approximate (delta-Eddington) 
code could be used to compute the upwelling spectral flux at the top of the 
atmosphere for the same surface - sun angle combinations but for a wide range 
of atmospheric conditions. Then by combining the fluxes with the appropriate 
angular directional models the complete spectral radiance models for all the 
desired atmosphere-surface-sun angle combinations may be approximated. This 
approach minimizes the use of the finite difference code which is very expensive 
in terms of computer costs. The final products therefore consist of a set of 
data files of spectral bidirectional models for each surface-sun angle com- 
bination and a set of files containing the spectral flux models. Together 
these files form a data base which may be used to provide the necessary ra- 
diances for the shortwave spectral correction algorithms. 

This report briefly summarizes the work which resulted in the above pro- 
ducts. It begins with a discussion of the development of a case list which 
was based on the availability of data (both surface reflectance and atmospheric 
data) and on a desire to minimize the number of full radiance calculations 



required due to the expense involved in running the finite difference code. 
This is followed by a brief description of the finite difference program, FD, 
with an emphasis on the implementation of the spectral, bidirectional reflect- 
ance data. Next is a discussion of the use of delta-Eddington code to compute 
the upwelling flux at TOA. The report ends with a description and examples of 
the data products which have been delivered for use in the scanner spectral 
correction algorithms. 

A follow on effort is planned in which all of the details of the computa- 
tions involved here will be presented and documented in a volume of the ERBE 
Science Reference Manual. 
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II. DETERMINATION OF SURFACE-ATMOSPHERE CASE STUDIES 


Regardless of the efficiency of the method used, any comprehensive computa 
tion of the radiance distribution at the top of the Earth’s atmosphere is a 
costly endeavor in terms of computer time and resources. The processes in- 
volved in describing the transfer of solar radiation through a variable atmo- 
sphere, the reflection by a wide variety of surfaces and the dependencies of 
these processes upon wavelength and illumination-viewing geometry make such a 
computation difficult at best. 

The requirement here is that just such a computation be made in order to 
provide a global data base of radiances for use in deriving the ERBE Scanner 
spectral correction coefficients. This data base needs to include examples of 
as many of the most common scene types visible from space as is practical. In 
addition, atmospheric variations can have important effects on the radiance dis 
tribution and so it is desirable to include as many typical atmospheres as pos- 
sible in the data base. Add to this the variation with solar zenith angle and 
it becomes evident that the number of cases can grow very quickly. With the 
large number of cases and the expense involved in making the radiance calcula- 
tions it becomes important to judiciously select the cases to be run. 

A. Initial Approach : 

This rationale provided the motivation for the initial study in this work, 
namely to develop a comprehensive set of case studies (i. e., scene type - atmo 
sphere combinations) and still keep the number of full radiance calculations 
to a minimum. The initial approach was to examine the spectra of the upwelling 
flux cases which were similar to the required data sets and to establish some 
measure of the similarity between differing scenes. For example, if it could 
be shown that a forest surface and a grass surface produced similar upwelling 



flux spectra under a range of atmospheric conditions, then, perhaps, the two 
scene types could be combined into one. If further similarities could be es- 
tablished in their angular responses, then an argument could be made to combine 
the two scene types into one for the radiance computations. 

An earlier study had produced upwelling flux spectra over the required 
spectral range, 0.25 - 5 pm for a number of surface atmosphere cases which were 
similar to those required in this study. The surfaces considered included ocean, 
forest, grass, dark soil, light soil and snow. Spectral albedo data for these 
surfaces was obtained from a report by Tarnopolskiy (Ref. 1). The atmospheric 
data came from MSL report #10 (Ref. 2) and include measurements of the vertical 
profiles of pressure, temperature, water vapor ans ozone over 106 meteorologi- 
cal stations throughout the world. The radiative transfer calculations were 
performed using a delta-Eddington approximation which was based upon the work 
of Shettle and Weinman (Ref. 3). Examples of these spectra are shown in Figure 1 

The criterion used for comparing these spectra was a "similarity factor" 
defined as follows: 

5 pm 

f M U) M U) d\ 

S i . = — y 

1 ’ J r 5pm „ 5pm 7 

f MU\) d\ f M Z U) d\ 

Uo J 0 J 

2 

where M^(\) is the spectral flux in W/m for the case designated by subscript 

i and \ is the wavelength. This quantity has a maximum value of 1.0 when 

M^(\) = Mj(\). Combinations of cases (i, j) which are similar therefore 

will have S. . z 1.0. 
i>J 
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This calculation was performed for a large number of the cases described 
above and results showed a high degree of similarity only for cases which had 
the same surface but with only slightly different atmospheric conditions. The 
decision was then made to use a different approach for minimizing the number 
of full radiance runs. 

B. Spectral Bidirectional Angular-Spectral Flux Approach : 

The method decided upon by which to form the required radiance models 
is based in the relationship between the radiance L and the flux (or radiant 
exitance) M through a bidirectional angular function R. This relationship 
may be defined as follows: 

R \ ( 0 O’ 9 > = 11 \ ( 0 O’ 05 /M x ( 0 o } (2) 

where (9 q, 0) is the spectral radiance in the viewing direction defined 

by zenith and azimuth angles 0 and <f>, respectively, resulting from illumina- 
tion by the sun at solar zenith angle 0^. Since the spectral flux M^(0^) is 
given by 


2t r it/ 2 

M.(0 n ) = / I L ( 0 n , 0, <j>) cos0 sin0 d0 d<J> , (3) 

\ u 0 0 * u 


a normalization condition may be written 


2tt tt/2 

tt J / R (0 n , 0, <j>) cos0 sin0 d0 d<J> = 1. (4) 

0 0 A 


Using the definitions given above, spectral bidirectional angular models 
may be derived from radiance simulations for the wide range of required scene 
types, but for a limited set of mean atmospheric conditions. Since the spectral 
flux is relatively inexpensive to compute (using a delta-Eddington model) a set 
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of spectral flux models could be developed which would embody a great number 
of atmospheric variations. Then by combining the bidirectional angular models 
with the spectral flux models the required radiances could be approximated. This 
approach limits the full radiance caculations to those for only mean atmospheric 
conditions and allows for the bulk of the atmospheric variations to be simula- 
ted by using the relatively inexpensive flux calculations. This was the approach 
adopted. 

C. Surface Reflectance Data ; 

The reflection properties of typical Earth surfaces are very complicated 
since they depend upon so many variables: wavelength, illumination - viewing 

geometry, time, etc. For this reason and because of the variety of surfaces 
and the problems involved with measuring these properties, sources of spectral, 
bidirectional reflection data are very sparse. One source of data was found 
to be very useful for the purposes required here, namely that of Kriebel 
(Ref. 4). Kriebel presents tables of spectral, bidirectional reflectance 
factors for four vegetated surfaces: savannah, bog, coniferous forest and 

pasture land, at seven narrow spectral regions (0.43, 0.52, 0.61, 0.87, 1.24, 

1.66 and 2.20 pm). The tables are presented for cones of reflection formed 
by dividing the hemisphere into 30° azimuth angle increments and 10° zenith 
angle increments. 

A second source of bidirectional reflectance data was a report by Davis 
and Cox (Ref. 5). In this report they present the results of an experiment 
in which measurements of the reflected radiance field over various atmospheric 
scene types were made using a multi-detector "bug-eye instrument" mounted on 
an aircraft platform and flown at an altitude of 10—12 km. The instrument 
had a broad-band spectral response ranging from 0.3 to 1.1 pm. The bidirec- 
tional angular models which they derived from their results were presented as 
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coefficients of a modified spherical harmonics expansion of the measured 
patterns. The data they presented for the desert sands of the Saudi Arabian 
Empty Quarters were used in this study. 

The above sources of data are limited in their spectral resolution. There- 
fore, in order to obtain reasonable representations of the wavelength variations 
of the surfaces of interest, it was decided to use some other sources of spec- 
tral albedo data in conjunction with the bidirectional angular data described 
above. Tarnopolskiy (Ref. 1) provides spectral albedo curves for some of the 
same surfaces for which the bidirectional data is available, namely, desert 
sands, coniferous forest and grass or pasture land. Spectral albedo data for 
snow and ice was obtained from a paper by Wiscombe & Warren (Ref. 6). By re- 
normalizing the spectral albedo data to the broad-band hemispherical albedo 
value computed from the bidirectional data, it was possible to approximate a 
spectral, bidirectional reflection model for several land surfaces: desert, 

forest and grass. Detailed spectral reflection data for the savannah surface 
was not available and so the coarse spectral resolution provided by Kriebel 
was used. In like manner, detailed angular reflection data for ice or snow 
surfaces was not available and so a Lambertian assumption was made for these 
cases. 

The ocean surface is somewhat different from typical land surfaces due to 
its pronounced specular reflection peak. This reflection of the directly in- 
cident beam may be described by a Fresnel law which is given in Siegel and 
Howell (Ref. 7). The diffuse reflection is assumed to be Lambertian and have 
a spectral variation as given by Morel (Ref. 8) for "blue (desert) water". 

Table 1 enumerates the surfaces decided upon for simulation in this 
study and lists the sources of the data used. 
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D. Atmospheric Data : 


As stated earlier, the full radiance calculations were to be done for 
"mean" atmospheric conditions, with simulations of the effects of detailed 
atmospheric variations being reserved for spectral flux calculations. The data 
for use in the radiance calculations comes from a report by McClatchey et al. 
(Ref. 9) and are for climatological conditions as follows: 1962 standard, 

tropical, midlatitude summer, midlatitude winter, subarctic (designated as 
"polar" for this work) summer and subarctic winter. The data consists of ver- 
tical distributions of the pressure, temperature, density, water vapor and 
ozone concentrations . 

The aerosol models used employ log-normal size distributions with para- 
meters and compositions appropriate to continental and maritime aerosols as 
specified in Deepak and Gerber (Ref. 10). The vertical distribution of aero- 
sols came from McClatchey et al. (Ref. 9) and is scaled to correspond to ground 
level visibilities of 5 km and 23 km. 

The data used to simulate atmospheric variability in the spectral flux 
models was developed from a data set which was obtained from the Tiros Opera- 
tional Vertical Sounder (TOVS) (Ref. 11). Data files were assembled by S. Gupta 
which provided monthly averages of the surface temperature, pressure, ozone 
burden, water vapor burden, cloud top height and vertical distribution of tem- 
perature for a 12 month period (July 1981 through June 1982) over 11 different 
sites. It seems reasonable to assume that these sites would have surface 
reflectance properties which could be approximated by the reflectance data 
described earlier. Table 2 enumerates these TOVS sites and the assumed scene 
type used to describe its reflectance properties in this work. 

The obscuration of the Earth* s surface by clouds provides another case of 
interest in providing a global data base of radiances. In the FD model clouds 
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are treated as a homogeneous aerosol layer with size distribution parameters 
as given in Stephens (Ref. 12). It was decided to simulate low, middle and 
high altitude clouds with optical depths of approximately 6, 12 and 20 respec- 
tively so that a range of effects could be examined. Clouds with tops below 
2.5 km are considered low; between 2.5 and 5.0 km, middle; and above 5.0 km, 
high. Stephen’s altostratus cloud droplet distribution parameters were used 
for the high and middle clouds and the stratocumulus II parameters used for 
the low clouds. The radiance field for these clouds over ocean and land (grass) 
surfaces were calculated as part of the data base. 

E. Case List : 

The resulting list of cases to be simulated using the FD radiance model 
are enumerated in Table 3. Four solar zenith angles were selected from the FD 
model’s angular grid for all but the polar winter scenes where the sun remains 
close to the horizon. Only the two largest solar zenith angles were used for 
these cases. The total number of radiance cases is 102. 

Table 4 lists the cases decided upon for spectral flux calculations. Both 
clear and cloudy cases were simulated except for the desert and ice/snow scenes. 
The cloud height for cloudy simulations was assigned according to the TOVS data. 
The solar zenith angles employed here were selected to match requirements of 
the spectral correction software. Considering each solar zenith angle - monthly 
average atmosphere - scene type combination as a unique case, the total number 
of spectral flux cases comes to 1080. 
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III. RADIATIVE TRANSFER MODELS 


As discussed above, the approach decided upon requires that full radiance 
calculations be done for a range of scene types and for mean, climatological 
atmospheric conditions so that spectral bidirectional models can be developed. 
These models, when combined with spectral flux models for the same scene types 
but for detailed atmospheric variations, could then be used to reconstruct a 
wide range of surface-atmosphere radiances. 

A. The Finite Difference Model : 

The radiative tranfer model used to compute the radiance field at the top 
of the atmosphere employs a finite difference technique to solve the radiative 
transfer equation for the Earth-atmosphere system. The model was developed by 
Barkstrom (Ref. 13) for solving the radiative transfer problem for a plane - 
parallel atmosphere with detailed vertical profiles of the physical and optical 
properties of the important atmospheric gases and aerosols. The model was 
enhanced by Suttles (Ref. 14) to include azimuthal variations, spectral varations 
for the entire solar wavelength rnage and specular in addition to diffuse re- 
flection at the surface. The model has been shown to yield results which are 
in good agreement with satelitte radiance observations (Ref. 14). 

The spectral range for which the program is set-up extends from 0.25 to 
5.0 ym. This wavelength region is broken up into 73 individual wavelength grid 
points at which the monochromatic radiative transfer equation is solved. This 
spectral region contains several strongly absorbing water vapor bands and to 
account for the absorption in these bands an exponential sum fit is used as 
described in Stephens (Ref. 12). This fit requires the solution to be repeated 
a number of times for each of the absorbing wavelengths which brings the total 
number of wavelengths in the spectral grid to 96. 
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The angular grid consists of 15 zenith angles and 11 azimuthal angles. The 
azimuthal angles are measured relative to the principal plane (containing the 
ray from the sun to the target area) with 4> = 0° corresponding to forward scat- 
ter and <p = 180° to backscatter. The zenith angle grid is chosen in such a man- 
ner as to minimize the error which results from the angular quadrature used in 
the solution. Therefore the angles used are unevenly spaced with regions of 
greatest variability in the phase function having more points. Table 5 gives 
a list of the zenith angles, the cosines of the zenith angles and the azimuthal 
angles in the solution grid. 

The depth grid is likewise determined by the solution technique and so the 
number of depth points depends upon the optical depth of the atmosphere and the 
accuracy of solution required. For the cases involved in this work the number 
of depth points was 49 for clear cases and 77 for clouds. Details on the selec- 
tion of angular and depth grids may be found in Barkstrom (Ref. 13) and Suttles 
(Ref. 14). 

Treatment of the azimuthal variation involves use of an expansion of the 
radiance field in an M-term Fourier cosine series in the azimuthal angle, 

M 

L (u, <f>) = l L® (y) cos m<|> 

K m=0 K 


In this formulation the lower boundary condition may be written 



where the first term represents the specular reflection at the surface, the 

second term represents the diffuse reflection of diffusely incident radiation 

and the third term represents the diffuse reflection of the directly incident 

solar beam. The factor p® (-y; y') is the Fourier cosine coefficient of the 

A 

spectral, bidirectional reflectance of the surface whose values are given by 
e 2 tt 

p” (y; U') = 27 / P x (u> ♦ » <t> r ) COS m( <f> - <J>') d(<J> - <j)') (7) 

where 

! 1 , if m = 0 

( 8 ) 

2, if m > 0 

and p (p, 4> ; p T , <f> ! ) are the diffuse spectral bidirectional reflectance data. 

A 

Evaluation of these coefficients was accomplished through the use of a 
fast Fourier transform routine available on the NASA LaRC NOS Mathematical 
Library (Ref. 15). The bidirectional reflectance data for the land surfaces 
described earlier (forest, grass and savannah from Kriebel and desert from 
Davis and Cox) were used to assemble a data base of spectral Fourier coeffic- 
ients which could then be used in the FD program. The data base was built by 
computing 128 terms in the Fourier, series for each of the elements in the 15 
by 15 angular matrix defined by p™ (p; p ! ). This was repeated for each wave- 

A 

length in the original data sets. Since the desert data of Davis and Cox was 
developed from a broad-band measurement, only one wavelength index was re- 
quired. The amount of data which resulted was rather large and therefore in 
order to facilitate access to the data by the FD program, the data base was 
stored on a random access file which allows for much more rapid retrieval than 
the standard sequential access storage. 
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The specular reflection term in equation (6) was assumed to be negligi- 
ble, except for the ocean surface in which it is a major contributor. The 
specular reflection from the ocean surface was assumed to obey a Fresnel re- 
flection law which is described in Siegel and Howell (Ref. 7). 

Figures 2 through 5 show some examples of the normalized spectral radi- 
ance reflected up at the top of the atmosphere as calculated using the FD 
model. Radiances are normalized to the value of the spectral radiant exitance 
at the top of the atmosphere. The radiance field is assumed to be symmetrical 
with respect to the principal plane. The spectrally integrated radiances 
(labeled as \ = TOTAL in the figures) are calculated numerically using 


"TOTAL 


■ L 1 

n=2 


+ L n-1> 


< X n- 


V-l> 


(9) 


where L^ is the monochromatic solution of the transfer equation for wave- 
length 

B. The Delta-Eddington Model : 

The delta-Eddington approximation is considered to be an efficient and 
reasonably accurate approach for radiative transfer calculations of the spectral 
flux. The method, as developed hy Shettle and Weinman (Ref. 3) and by Joseph 
et al. (Ref. 16), is based on Eddington's approximation in which the azimuthal 
average radiance is written as a two term expansion in the cosine of the zenith 
angle. 


L(t, y) = L q (t) + y L^t). 


( 10 ) 


The next step is to assume that the phase function for scattering may be 
approximated by using a Dirac delta function for the forward scatter peak and 
a two term expansion in the scattering angle a: 
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P(cosa) = 2 f 6(1 - cosa) + (1 - f) (1 + 3g' cosa) 


(ID 


where f is the fraction scattered into the forward direction and g' is the 
mean value of the cosine of the scattering angle (or the asymmetry factor). 
Under these assumptions and approximations, the radiative transfer equation 
for an atmosphere consisting of homogeneous layers, can be reduced to a pair 
of coupled ordinary differential equations which can then be solved in closed 
form. 

Software to perform these calculations had already been set up by Suttles 
and Tiwari (Ref. 17) with subroutines detailing the solar spectral irradiance 
and the Earth-atmosphere properties shared in common with the FD model. This 
gave further justification to its use in computing the spectral fluxes for 
radiance data base. 
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IV. DATA PRODUCTS 


The primary goal of this research effort is to provide a data base of 
short-wave, spectral radiance models which are to be used in algorithms to 
correct for spectral effects caused by the ERBE Scanner Instrument. As dis- 
cussed above, it was decided to use a detailed finite difference model to de- 
rive a set of spectral bidirectional angular models for the required Earth 
scene types and to combine these with a set of spectral flux models computed 
using a del ta-Eddington model to simulate typical atmospheric variabilities. 
These data files were prepared and delivered to the ERBE Inversion Working 
Group in October, 1984. 

A. Spectral Bidirectional Angular Models : 

The spectral correction algorithms required that the radiance data base 
be stored using the angular bin structure shown in Table 6. FD grid points 
located near the centers of the solar zenith angle bins were considered to be 
representative for the entire bin. The same was true for the 5 relative 
azimuth angle bins, i. e., the value for the angle grid point at the center 
was considered to be representative of the values over the entire bin. To 
match the viewing zenith angle bins however, it was decided to combine the 
radiance values at adjacent zenith angle grid points and to compute an average 
radiance value for the bin. Weighted sums of the radiances across the FD angle 
ranges as indicated in Table 6 were computed with the resulting average value 
being stored as the values for the appropriate bins. The radiance values were 
all normalized to the values of the spectral flux as indicated in equation (2) 
to obtain the bidirectional angular models. 
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These resulting models were stored on data files identified by the scene 
type and solar and solar zenith angle. Bidirectional angular model values 
were stored as 4 x 5 arrays (4 viewing zenith bins by 5 azimuth bins) for each 
of the 73 wavelength values from 0.25 to 5.0 pm. An example is shown in Figure 
6 . 

B. Spectral Flux Models : 

The spectral flux models were assembled in order to provide variability in 
the radiance models due to atmospheric variations. With atmospheric data com- 
ing from the TOVS models as discussed earlier, the spectral flux data files 
comprised 73 spectral flux values for 4 solar zenith angles (only 2 solar zenith 
angles for polar winter cases) and for 12 months of atmospheric data. Figure 7 
shows an example of one month of data stored a one of the 29 data files which 
makeup the spectral flux models. 
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V. CONCLUDING REMARKS 


The models of spectral radiance described in this report have been created 
with a specific purpose in mind, namely to aid in the determination of the ERBE 
Scanner spectral correction factors. The data base created herein provides 
the capability to simulate radiance fields over the most common Earth scene 
types for a large number of atmospheric conditions. The atmospheric properties 
have been derived from recent measurements and the bidirectional reflectance 
properties employed are as realistic as are currently available. These models 
may also form the basis for a number of studies involving the angular effects 
of surface reflection and effects of atmospheric variability. 

The next level of detail which would be desired in further studies would 
involve the inclusion of wind-roughened ocean surface effects and the inclusion 
of embedded clouds within the atmosphere. These are being planned for the near 
future. 
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TABLE 1. SURFACE TYPES AND REFLECTION DATA SOURCES 


SURFACE 

SPECTRAL ALBEDO DATA 

ANGULAR BIDIRECTIONAL MODEL 

Forest 

Tarnopolskiy (Ref. 1) 
(forest winter species) 

Kriebel (Ref. 4) 
(coniferous forest) 

Grass 

Tarnopolskiy 

(dry valley meadows) 

Kriebel 

(pasture land) 

Savannah 

Kriebel 

Kribel 

Desert 

Tarnopolskiy 

(sand, desert outcrop, rocks) 

Davis & Cox (Ref. 5) 

(Saudi Arabian Empty Quarter) 

Ice/Snow 

Wiscombe & Warren (Ref. 6) 
(dry, old snow) 

(assumed Lambertian) 

Ocean 

Morel (Ref. 8) 

(blue, desert water) 

Siegel & Howell (Ref. 7) 
(Fresnel Law) 
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TABLE 2. TOVS METEOROLOGICAL DATA SITES AND ASSUMED 
SCENE TYPES FOR SIMULATION 


. SET NO. 

SITE DESCRIPTION 

LAT - LON 

SCENE TYPE(s) 

TO VI 

Tropical Ocean - Pacific 

20° N - 130° W 

Ocean 

TOV2 

Tropical Ocean - Indian 

0'lat - 90 ° E 

Ocean 

TOV3 

Tropical Land - Amazon Jungle 

5° S - 60° W 

Forest, Grass 

TOV4 

Tropical Desert - Sahara 

22.5°N - 25° E 

Desert, Savannah 

TOV5 

Tropical Desert - Saudi Arabia 

20° N - 50° E 

Desert 

TOV6 

Midlat Ocean - N. Atlantic 

45° N - 40° W 

Ocean 

TOV7 

Midlat Land - Midwest U. S. 

40° N - 95° W 

Grass, Savannah 

TOV8 

Midlat Snow - W. Canada 

50° N - 110° W 

Snow (Winter, Spring) 
Forest (Summer, Fall) 
Grass (Summer, Fall) 

TOV9 

Polar Ocean - North Sea 

70° N - 0° 

/u LON 

Ocean 

TOVIO 

Polar Land - Siberia 

67° N - 100° E 

Snow (Winter, Spring) 
Forest (Summer, Fall) 

TOV11 

Polar Ice/Snow - Greenland 

75° N - 40° W 

Ice/Snow 











TABLE 3. RADIANCE CASE LIST 


SCENE 

ATMOSPHERE 

SOLAR ZENITH ANGLES (DEGREES) 

Ocean 

Tropical 

— Mari , 

V 

= 

5 

33.4, 

55.6, 

69.7, 

84.0 


Midlat 

Summer fl 

V 

= 

5 

II 

II 

II 

II 


Midlat 

Winter " 

V 

= 

23 

II 

II 

It 

II 


Polar 

Summer " 

V 

= 

23 

It 

II 

II 

II 


Polar 

Winter " 

V 

= 

23 

- 


69.7, 

84.0 

Grass 

Tropical 

Cont, 

V 

= 

5 

33.4, 

55.6, 

69.7, 

84.0 


Midlat 

Summer " 

V 

= 

5 

II 

II 

II 

II 


Midlat 

Winter " 

V 

= 

23 

II 

It 

II 

II 

Forest 

Tropical 

Cont , 

V 

= 

5 

33.4, 

55.6, 

69.7, 

84.0 


Midlat 

Summer 11 

V 

= 

5 

II 

II 

II 

II 


Midlat 

Winter !f 

V 

= 

23 

II 

II 

II 

It 


Polar 

Summer 11 

V 

= 

23 

II 

II 

II 

II 


Polar 

Winter " 

V 

= 

23 

- 

- 

69.7, 

84.0 

Savannah 

Tropical 

— Cont , 

V 

s 

5 

33.4, 

55.6, 

69.7, 

84.0 


Midlat 

Summer 11 

V 

= 

5 

It 

It 

II 

II 


Midlat 

Winter " 

V 

= 

23 

II 

II 

II 

It 

Desert 

Tropical 

— Cont , 

V 

= 

23 

33.4, 

55.6, 

69.7, 

84.0 

Ice/Snow 

Midlat 

Summer Cont , 

V 

— 

5 

33.4, 

55.6, 

69.7, 

84.0 


Midlat 

Winter " 

V 

= 

23 

II 

II 

II 

II 


Polar 

Summer " 

V 

= 

23 

It 

II 

II 

II 


Polar 

Winter " 

V 

= 

23 

— 


II 

II 

Cloud/Ocean 

Low Cloud , > x = 6 




33.4, 

55.6, 

69.7, 

84.0 


Middle Cloud, x= 12 




II 

II 

II 

II 


High Cloud, x= 20 




II 

II 

It 

II 

Cloud/Land 

Low Cloud, t = 6 




33.4, 

55.6, 

69.7, 

84.0 


Middle Cloud, T= 12 




It 

II 

II 

II 


High Cloud, t= 20 




II 

II 

II 

II 
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TABLE 4. SPECTRAL FLUX CASE LIST 


SCENE 

LAT. ZONE 

MONTHS 

TOVS DATA 3 

SOLAR ZENITH ANGLES (DEGREES) 

Ocean 3 

Tropical 

Jul 

_ 

Jun 


29.0, 

51.3, 

68.4, 

82.8 


Tropical 

Jul 

- 

Jun 

l 

29.0, 

51.3, 

68.4, 

82.8 


Midi at 

Jul 

- 

Jun 


II 

If 

tf 

o 

If 

o 


Polar 

Jul 


Jun 

T0V9 

- 

- 

68.4 , 

82.8 

Grass 3 

Tropical 





29.0, 

51.3, 

68.4, 

82.8 


Midlat 





ft 

ft 

If 

tf 


Midi at 




TOV8 

It 

ft 

If 

If 

Forest 3 

Tropical 

Jul 

_ 

Jun 


29.0, 

51.3, 

68.4, 

82.8 


Midlat 

Apr 

- 

Sep 


II 

It 

If 

fl 


Polar 

Apr 

- 

Sep 

TOVIO 

It 

tf 

ft 

It 

Savannah 3 

Tropical 

Jul 

_ 

Jun 

T0V4 

29.0, 

51.3, 

68.4, 

82.8 


Midlat 

Apr 

- 

Sep 

T0V7 

It 

ft 

ft 

II 

Desert* 3 

Tropical 

Jul 

_ 

Jun 


29.0, 

51.3, 

68.4, 

82.8 


Tropical 

Jul 


Jun 

■ 

II 

tf 

If 

tl 

Ice/Snow^ 

Midlat 

Oct 

— 

Mar 


29.0, 

51.3, 

68.4, 

82.8 


Polar 

Oct 

- 

Mar 



- 

68.4, 

82.8 


Polar 

Apr 

“ 

Sep 

T0V11 

29.0, 

51.3, 

68.4, 

82.8 


Polar 

Oct 


Mar 

T0V11 

— 

— 

68.4, 

82.8 


Both clear and cloudy cases 
^Clear cases only 
C See Table 2 
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TABLE 5. ANGULAR SOLUTION GRID FOR FD MODEL 


ZENITH ANGLES, 0 (DEGREES) 

COS (0) 

88.790 

0.0211 

84.039 

0.1039 

76.997 

0.2250 

69.748 

0.3461 

64.603 

0.4289 

61.741 

0.4735 

55.571 

0.5654 

45.573 

0.7000 

33.424 

0.8346 

22.098 

0.9265 

17.471 

0.9539 

14.592 

0.9677 

10.226 

0.9841 

5.541 

0.9953 

1.793 

0.9995 


AZIMUTHAL ANGLES, <j> (DEGREES) 


0.0 

7.5 

30.0 

37.5 

60.0 

90.0 

120.0 

142.5 

150.0 

172.5 

180.0 





TABLE 6. ANGULAR BINS FOR SPECTRAL BIDIRECTIONAL MODELS 


Solar Zenith 
Angle Bins 


y 

0 

0 q (degrees) 




From 

To 

From 

To 

Diff . 

Center 

Nearest 
FD Point 

1.00 

0.75 

0 

41*4 

41.4 

29.0 

33.4 

0.75 

0.50 

41.4 

60.0 

18.6 

51.3 

55.6 

0.50 

0.25 

60.0 

75.5 

15.5 

68.4 

69.7 

0.25 

0.00 

75.5 

90.0 

14.5 

82.8 

84.0 



0 (degrees) 


From 

To 

Diff. 

Center 

FD Angle 
Range 

Indices 


0 

30 

30 

15.0 

1.8 - 22.4 

15 - 10 

Viewing Zenith 

30 

45 

15 

37.5 

33.4 

9 

Angle Bins 

45 

60 

15 

52.5 

45.6 - 55.6 

8-7 


60 

90 

30 

75.0 

61.7 - 88.8 

6 - 1 



<f> (degrees) 



From 

To 

Diff. 

Center 


0 

15 

15 

7.5 

Relative Azimuth 

15 

60 

45 

37.5 

Angle Bins 

60 

120 

60 

90.0 


120 

165 

45 

142.5 


165 

180 

15 

172.5 
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Figure 1. Sample of Top of Atmosphere Flux Spectra Used in Similarity Study. 
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0MIDW33 — OCEAN, MIDLAT WINTER, MARITIME AERO, V=23., SUN=33 DEG 

WAVELENGTH = .55 

C0NT0UR FR0M .70000 T0 20.000 C0NT0UR 1NTERVRL 0F SPECIRL TENSION 0F 3.0000 


90.00 


M. =131.10 W/m /um 


ZENITH ANGLE, 9 

Figure 2. Normalized Top of Atmosphere Radiances for Ocean Surface, Midlatitude Winter, 0 


(a) \ = 0.55 pm 






MARITIME AERO, V=23„ SUN=33 DEG. 
= 1.00 


:rvrl bf specirl tens i bn bf 3.0000 

2 

M = 29.34 W/m /pm 



ean Surface, Midlatitude Winter, 9 = 33.4°. 

o 

pm 


AZIMUTH ANGLE 


ZZ — OCEAN, 

MIDLAT 

WINTER, MARITIME AERO, V= 

=23 

SUN=33 DEG 


WAVELENGTH = TOTAL 



C0NT0UR FR0H .70000 

T0 20.000 

C0NT0UR INTERVAL 0F SPECIAL TENSION 

0F 

3.0000 



90.00 

M = 

92.78 W/m 2 


h/WA 


55.00 


'S.uu 03. UU 33. LU 13. uu 15.00 35.00 55.00 

ZENITH ANGLE, 0 

Figure 2. Normalized Top of Atmosphere Radiances for Ocean Surface, Midlatitude Winter 

(c) \ = TOTAL 




FMIDW33 — FOREST, MIDLAT WINTER, CONT AERO, V=23„ SUN=33 DEG 

WAVELENGTH = 1.00 

C0NT0UR FR0M .70000 T0 20.000 C0NT0UR INTERVAL 0F SPECIAL TENSION 0F 3.0000 


90.00 


M. =? 111.94 W/m /lim 


ZENITH ANGLE, 8 


Figure 3. Same as Figure 2, Except for Forest Surface. 


35.00 55.00 


(a) A = 0.55 ym 





FMIDW33 — 

FOREST, 

, MIDLAT WINTER, 

CONT 

AERO, 

V=23-, 

SUN =33 

DEG 



WAVELENGTH 

= .55 





C0NT0UR FR0M 

.70000 

T0 20.000 C0NT0UR INTERVAL 0F 

SPECIAL 

TENSI0N 0F 

3.0000 




90.00 




M. = 76.99 

2 

W/m /ym 


ZENITH ANGLE, 0 

Figure 3. Same as Figure 2, Except for Forest Surface. 

(b) A = 1.00 um 



FMIDW33 — 

FOREST, 

, MIDLAT WINTER, CONT AERO, 

V=23., SUN=33 DEG 



WAVELENGTH = TOTAL 


CONTOUR FROM 

.70000 

T0 20.000 CONTOUR INTERVAL OF SPECIAL 

TENSION OF 3.0000 



90.00 

M = 109.04 W/m 2 




15.00 15.1 



ZENITH ANGLE, 0 

Figure 3- Same as Figure 2, 

Except 

for Forest Surface. 



(c) X = TOTAL 


35.00 






IMIDW33 — ICE-SNOW, MIDLAT WINTER, CONT AERO, V=23., SUN=33 DEG 

WAVELENGTH = .55 

C0NT0UR FR0M .70000 T0 20.000 C0NT0UR INTERVAL 0F SPECIAL TENSION 0F 3.0000 


90.00 


= 107.49 W/m /]_lm 



75.00 55.00 35.00 


15.00 


15.00 35.00 55.00 75.00 


ZENITH ANGLE, 6 


Figure 4. Same as Figure 2, Except for Ice/Snow Surface. 


(a) X - 0.55 ym 








IMIDW33 — ICE-SNOW, MIDLAT WINTER, CONT AERO, V=23„ SUN=33 DEG 

WAVELENGTH = TOTAL 

C0NT0UR FR0M .70000 T0 20.000 C0NT0UR INTER VRL 0F SREC1RL TENSION 0F 3.0000 


90.00 


M = 619.92 W/m 



75.00 


55.00 


35.00 


15.00 


15.00 


35.00 


55.00 


ZENITH ANGLE, 0 


Figure 4. Same as Figure 2, Except for Ice/Snow Surface. 
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(c) A = TOTAL 



J72.S 


DTRP33C 


DESERT, TROPICAL ATM, CONT AERO, V=23., SUN=33 DEG 


WAVELENGTH = .55 

C0NT0UR FROM .70000 T0 20.000 CONTOUR INTERVAL OF SPECIRL TENSION OF 3.0000 



75.00 


55.00 


35.00 


15.00 15.00 

ZENITH ANGLE, 0 


35.00 


55.00 


75.00 


Figure 5. Same as Figure 2, Except for Desert Surface, Tropic. 


(a) X = 0.55 ym 


AZIMUTH ANGLE 



* 72.5 


* 1 

DTRP33C DESERT, TROPICAL ATM, CONT AERO, V=23., SUN=33 DEG 

WAVELENGTH =1.00 

C0NT0UR Ff?0M .70000 T0 20.000 C0NT0UR INTERVAL 0F SPECIAL TENSION 0F 3.0000 


90.00 M x = 217.99 W/m 2 /ym 



(b) X = 1.00 y m 


AZIMUTH ANGLE 






SPECTRAL BIDIRECTIONAL MODEL, CASE: 

FTRQP33 FOREST, TROPICAL ATM, CONT AERO, V=5. , SUN=33 


WAVELENGTH = .250 MICRONS 


ZENITH 
AZIMUTH 
- 7.5 

37.5 
90. 0 
* 1 42. 5 
172.5 

> 15.0 

. 6861 10E+00 
. 700258E+00 
. 763167E+00 
. 842676E+00 
. 8651 86E+00 

37 . 5 

. 665506E+00 
. 687770E+00 
. 781843E+00 
. 933609E+00 
. 980580E+00 

52.5 

. 797 1 84E+00 
. 792099E+00 
. 86 1 044E+00 
. 107174E+01 
. 1 1 4651E+01 

75.0 

. 149019E+01 
. 140367E+01 
. 13001 IE +01 
. 158257E+01 
. 171328E+01 

WAVELENGTH = . 260 

ZENITH > 15.0 

MICRONS 

37.5 

52.5 

75. 0 

AZIMUTH 

7. 5 

. 684320E+00 

. 659667E+00 

. 791936E+00 

. 1 49036E+01 

37.5 

. 699029E+00 

. 684989E+00 

. 787819E+00 

. 140411E+01 

90. 0 

. 763099E+00 

. 781236E+00 

. 859657E+00 

. 130244E+01 

142.5 

. 843674E+00 

. 934992E+00 

. 1 07302E+01 

. 158777E+01 

172.5 

. 86645 1E+00 

. 982453E+00 

. 1 1 4844E+01 

. 171936E+01 

WAVELENGTH = . 270 

ZENITH > 15.0 

MICRONS 

37.5 

52. 5 

75.0 

AZIMUTH 

7.5 

. 682094E+00 

. 622 1 05E+00 

. 726722E+00 

. 145346E+01 

37.5 

. 700439E+00 

. 664009E+00 

■ 735072E+00 

. 13756 IE +01 

90. 0 

. 777010E+00 

. 7B5093E+00 

. 844816E+00 

• 1 3 U 8 7 1 E + 0 1 

142.5 

. S70189E+00 

. 963337E+00 

. 109296E+01 

. 163243E+01 

172. 5 

. 896163E+00 

. 101693E+01 

. 1 17667E+01 

. 1 77420E+01 




WAVELENGTH = 4.500 

MICRONS 



ZENITH 

> 15.0 

37.5 , 

52. 5 

75.0 

AZIMUTH 





7.5 

. 668334E+00 

. 720128E+00 

. 887937E+00 

. 125278E+01 

37.5 

. 739086E+00 

. 906169E+00 

. 790538E+00 

. 105016E+01 

90 . 0 

. 767395E+00 

. 810163E+00 

. 8201 17E+00 

. 940876E+00 

142.5 

. 102191E+01 

. 127950E+01 

. 1 25204E+0 1 

. 1 10821E+01 

172. 5 

. 124898E+01 

. 1 90524E+01 

.. 1 96388E+0 1 

. 146403E+01 

WAVELENGTH = 5.000 

MICRONS 



ZENITH 

> 15.0 

37.5 

52 . 5 

75.0 

AZIMUTH 





7.5 

• 667563E+00 

. 721474E+00 

. 887580E+00 

. 121850E+01 

37.5 

. 740542E+00 

.8101 16E+00 

. 737862E+00 

. 10197 1E+01 

90 . 0 

. 768295E+00 

. 812877E+00 

. 820130E+00 

. 929309E+00 

142. 5 

. 102798E+01 

. 1 29201E+01 

. 1 26 1 57E+01 

. 1 1 0538E+0 1 

172. 5 

. 126156E+01 

. 193640E+01 

. 199613E+01 

. 1 47445E+0 1 


Figure 6. Sample of Spectral Bidirectional Angular Models. 
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JUL 


FOREST CONIFER SPECIES 


TROPIC 


V 


KM 


SPECTRAL FLUX UP AT TOA, W M-2 MICRON-1 


29.0 


SOLAR ZENITH ANGLE, DEGREES 


51.3 


68. 4 


82.8 


WAVELENGTH, 

MICRONS 

. 2500 OE +00 
. 26000 E +00 
. 27 000 E+ 00 
. 28000E+00 
. 29000E+00 
. 29500E+00 
. 30000E+00 
. 30500E+00 
. 31000E+00 
. 31500E+00 
. 32000E+00 
. 32500E+00 
. 33000E+00 
. 34000E+00 
. 35000E+00 
. 36000E+00 
. 37000E+00 
. 38000E+00 


. 13936E-01 
. 22489E-01 
. 47910E-01 
. 7436 IE— 01 
. 351 74E+00 
.71 737E+00 
. 1 0394E+0 1 
. 48334E+01 
. 27970E+02 
. 81541E+02 
. 13615E+03 
. 20635E+03 
. 25174E+03 
. 27 1 1 3E+03 
. 26607E+03 
. 24457E+03 
. 25390E+03 
. 2261 1E+03 


. 1 2068E-0 1 
. 194S3E-01 
. 41618E-01 
. 64936E— 01 
. 307B3E+00 
. 62731 E+00 
. 88944E+00 
. 31855E+01 
. 1 8826E+02 
. 61747E+02 
. 1 1023E+03 
. 17482E+03 
. 21913E+03 
. 24337E+03 
. 24297E+03 
. 22632E+03 
. 23792E+03 
. 21445E+03 


. 90663E— 02 
. 14644E-01 
. 31397E-01 
. 49485E-01 
. 23593E+00 
.48031 E+00 
. 67 1 82E+00 
. 17204E+01 
. 81870E+01 
. 32639E+02 
. 66B68E+02 
. 1 1673E+03 
. 15462E+03 
. 1 8 1 79E+03 
. 1 8628E+03 
. 1 7653E+03 
. 18851E+03 
. 17246E+03 


.4161 OE— 02 
. 67232E-02 
. 14452E— 01 
. 23157E-01 
. 11 381 E+00 
. 23220E+00 
. 32292E+00 
. 7 1 289E+00 
. 16312E+01 
. 52107E+01 
. 1 3956E+02 
. 32703E+02 
. 52108E+02 
. 73291 E+02 
. 79565E+02 
. 77150E+02 
. 83925E+02 
.78081 E+02 






. 18/00E+01 
. 20700E+0 1 
. 23000E+0 1 
. 27000E+0 1 
. 30000E+01 
. 31 000E+01 
. 32000E+0 1 
. 35000E+0 1 
. 40000E+01 
. 4 500 OE +01 
. 50000 E +01 

TOTAL FLUX 
ALBEDO 


. 26360E+01 
. 65643E+01 
. 51328E+01 
. 35959E— 01 
. 14108E+01 
. 12035E+01 
. 32358E— 01 
. 69290E+00 
. 43686E+00 
. 266 1 2E+00 
. 17041 E+00 

. 1 2688E+03 
. 1 1 1818 


. 21465E+01 
. 59432E+01 
. 45730E+01 
. 27 1 07E-0 1 
. 12618E+01 
. 10812E+01 
. 23845E— 01 
.62781 E+00 
. 39274E+00 
. 23785E+00 
. 15153E+00 

. L2437E+03 
. 153326 


. 1 3640E+0 1 
. 48215E+01 
. 36181E+01 
. 1 6220E-01 
. 1 0062E+01 
. 86913E+00 
. 1 3442E-01 
. 51218E+00 
. 3 1 659E+00 
. 18985E+00 
. 1 1994E+00 

. 10591E+03 
. 221757 


■ 37203E+UO 
. 27076E+01 
. 19200E+01 
. 70253E— 02 
. 54984E+00 
. 48666E+00 
. 48394E— 02 
. 29862E+00 
. 17995E+00 
. 10513E+00 
. 65040E— 01 

. 52683E+02 
. 323994 


Figure 7. Sample of Spectral Flux Models. 
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